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The local structural distortions in polydomain ferroelectric PbTiOs/SrTiOs superlattices are in- 
vestigated by means of high spatial and energy resolution electron energy loss spectroscopy combined 
with high angle annular dark field imaging. Local structural variations across the interfaces have 
been identified with unit cell resolution through the analysis of the energy loss near edge structure of 
the Ti-L2,3 and 0-K edges. Ab-imtio and multiplet calculations of the Ti-L2,3 edges provide unam- 
biguous evidence for an inhomogeneous polarization profile associated with the observed structural 
distortions across the superlattice. 



Complex oxide heterostructures offer a vast play- 
ground for exploring and combining the many functional 
properties of these fascinating materials arising from the 
subtle interplay between their charge, spin, orbital and 
lattice degrees of freedom. Bilayers, multilayers and su- 
perlattices composed of ultrathin oxide layers not only 
shed new light on our fundamental understanding of the 
constituent materials, but frequently reveal unexpected 
new phases at their interfacefP^. Superlattices composed 
of ferroelectric and paraelectric oxides have been the sub- 
ject of numerous studies, motivated by fundamental ques- 
tions about ferroelectric size effects, by the possibilities 
these artificially layered materials offer for tailoring their 
functional properties, and by the fascinating interface 
physics they display. Ultrafine period superlattices, com- 
posed of ferroelectric PbTiOa (PTO) and paraelectric 
SrTiOs (STO), for example, have been shown to exhibit 
improper ferroelectricity driven by the coupling of the 
polar and non-polar distortions at the interface^. More 
recently, regular ferroelectric nanodomains have been ob- 
served in PTO/STO superlattices with larger periodici- 
ties and were shown to be responsible for large enhance- 
ments in the effective dielectric constanlP. Such domains 
are expected to give rise to complex inhomogeneous 
structural distortions and polarization profilePEI, depart- 
ing from uniform polarization models, frequently used to 
describe the properties of ferroclectric/paraelectric su- 
perlattices in the absence of domains^. 

A microscopic insight into the local structure is thus 
key to understanding the behavior of these exiting arti- 
ficially layered materials, and here transmission electron 
microscopy (TEM), with the recent advances in spatial 
resolution, provides an invaluable tool. Individual ionic 
displacements within a single perovskite unit cell can now 
be identified from phase contrast images. The spatial 
resolution is high enough to determine the direction and 
even the magnitude of the local dipole moments in ferro- 
electric materialJ^^, and has recently enabled the direct 
verification of the existence of polarization rotation at 
domain walls in Pb(Zr,Ti)03 ferroelectric thin filmpl. 




FIG. 1. (a) HAADF-STEM image of the [010] projection of a 
(6|6) PTO/STO superlattice (b) X-ray reciprocal space scans 
around the (002) superlattice reflection show broad satellite 
peaks due to 180 ° ferroelectric domains. 



In this article, we focus on an alternative spectroscopic 
technique to study local ferroelectric distortions at the 
single unit-cell scale. Atomically-resolved electron energy 
loss spectra (EELS) and high angle annular dark field 
(HAADF) images have been obtained to access the local 
structur e!^ ^tf^ of polydomain PTO/STO superlattices. 
We show that the energy loss near edge fine structure 
(ELNES) of the Ti-L2^3-edge is highly sensitive to very 
small atomic displacements (< 0.1 A). Changes in the 
tetragonal distortion of the perovskite unit cell as small 
as 1% have been detected, revealing a variation of tetrag- 
onality within individual PTO layers and thus provid- 
ing direct evidence of local structural inhomogeneities in 
polydomain ferroelectric/paraelectric superlattices. Ab- 
initio and charge transfer multiplet calculations show dis- 
tinct differences between the spectra expected for polar 
and non-polar structures, demonstrating that the EELS 
technique has sufficient energy and spatial resolution to 
probe ferroelectricity at the unit-cell scale. 

We have studied a (6|6)2i superlattice with 21 bilay- 
ers of 6 unit cells of PTO and 6 unit cells of STO on 
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FIG. 2. (a) Ti-L2,3 and (b) O-K edges for STO substrate 
(black-upper) and 12 nm thick PTO layer (red-lower), (c) 
Schematic projection along the [010] zone-axis of the STO 
and PTO unit cells. The displacement of the atomic position 
is displayed as 5{Ti), 5{0l) and 5(02). 



top of a single crystalline STO substrate. Epitaxial top 
and bottom SrRuOa electrodes were also deposited in 
situ. Details of sample growth and electrical charac- 
terization can be found in Ref. |31 Atomically-resolved 
HAADF images and EELS spectra were acquired with 
an aberration-corrected scanning transmission electron 
microscope (STEM) Nion UltraSTEM 100^^. Figure la 
displays a HAADF image of the superlattice projected 
onto the (010) plane. The variation in HAADF signal 
intensity for each atomic column reflects the difference 
in atomic numbers between Pb (Z=82) and Sr (Z=38) 
cations^'', demonstrating the atomic sharpness of the in- 
terfaces. X-ray diffraction measurements reveal a peri- 
odic in-plane modulation attributed to 180° ferroelectric 
domains (Fig. lb). The domain satellites are observed in- 
dependent of the in-plane orientation of the sample, in- 
dicating an almost isotropic distribution of domain wall 
alignments. The domain periodicity ranges from 55 A 
(for domain walls along (100)) to 65 A (along (110)). 
This means that in general, the domain wall orientation 
will be random with respect to the imaging plane and the 
HAADF and EELS images will average over the domains 
within the thickness of the TEM specimen. 

As a reference for the characterization of the (6|6) su- 
perlattice, Ti-L2,3 and O-K EELS spectra were acquired 
from a 12 nm thick PTO layer and from the STO sub- 
strate (7 nm away from the interface). Figure 2a shows 
the Ti-L2,3 ELNES of the reference spectra for PTO and 
STO to be compared with those of the superlattice. The 
energy resolution (0.4 eV) achieved allows us to iden- 
tify slight variations between the STO and PTO Ti-L2_3 
ELNES that can be correlated with structural differences 
between the two oxides. While the unit cell of the STO 
substrate has cubic symmetry (a=6=c=3.905 A), the hy- 
bridization between the Pb 6s and the O 2p states in- 
duces a large tetragonal strairl^^ in the PTO unit cell 
(c/a=1.063) (Fig. 2c). Due to the reduction in symmetry 
from the cubic perovskite structure, the crystal field split- 
ting (CFS) of the Ti 3d e^ and t2g orbitals is expected 
to decrease, leading to the observed reduction in the La- 



edge energy splitting ( ALg)?^ from 2.3 eV in the STO 
substrate to 1.65 eV for the PTO layer (Fig.2a). The 
site-symmetry of the Ti atom in the TiOg octahedra is 
also reduced in the PTO unit cell, the corresponding Ti 
Lg-Cg peak showing an asymmetric broadening towards 
higher energies (see Fig. 2a), that has been attributed to 
the non-centrosymmetric location of the Ti cation in the 
TiOe octahedrai5 2i. 

Comparing the O-K ELNES of the STO substrate and 
the 12 nm thick PTO layer in Fig. 2b, clear differences 
between both spectra can be identified. The most pro- 
nounced one is the large shift observed in the position of 
the second peak (marked with an asterisk). This peak 
has been attributed to the hybridization of the O 2p 
with Sr 4d and with Pb 6sp states for STO and PTO, 
respectiveljl22H24|^ making it possible to distinguish be- 
tween Sr and Pb containing cells by identifying the po- 
sition of this peak. The analysis of the Ti-L2,3 and O-K 
edges therefore offers an excellent method for determin- 
ing both structural and chemical variations within the 
PTO/STO superlattices. 

Figure 3a shows a series of 512 individual EELS spec- 
tra acquired in spectrum- line mode with dwell time=125 
ms, displaying the Ti-L2,3 and the O-K edges across four 
interfaces of the (6|6) superlattice. The high spatial res- 
olution of STEM allows us to identify the core loss edges 
from individual atomic columns, while the dark field pro- 
file is simultaneously recorded along the scanning line 
(black line in Fig.3b). The Ti-L2,3 and O-K EELS spec- 
tra corresponding to 7 consecutive unit cells across one of 
the PTO-STO interfaces (numbered from 6 to 12 in the 
image) are displayed in Fig. 3d and 3e, respectively. Each 
spectrum is the sum of 22 individual spectra over one 
unit cell. In order to improve the energy resolution, the 
Ti-L2,3 spectra in Fig.3f have been deconvoluted using 
the Richardson-Lucy (R-L) approaclP^. 

Consistent with their large tetragonality, the Ti-L3 
splittings in the PTO layers are smaller than those in 
STO layers. The value for each unit cell along the super- 
lattice is displayed in Figure 3c. For better accuracy, the 
positions of the t2g and e^ peak maxima were obtained 
from the first derivatives of the EELS spectra. The split- 
ting value in the center of the PTO layers reaches 1.95 eV, 
significantly higher than the 1.65 eV obtained for the ref- 
erence 12 nm PTO layer, and increasing further on ap- 
proaching the interface. Interestingly, a gradual broad- 
ening of the L3-eg band is also observed when moving 
from the interface (spectrum 9) to the center of the PTO 
layer (spectrum 12). The spectral splitting of the STO 
layers is closer to the value of the STO substrate (2.3 
eV), but always around 50 — 100 meV smaller. In order 
to rule out the interdiffusion of Sr and Pb as the cause 
of the gradual variation in PTO, chemical profiles across 
the interface were obtained by multiple least squares fit- 
ting of the O-K edge signal from 536 to 540.8 eV with 
the two reference components (red and green spectra in 
Fig.3e). The resulting Sr and Pb distributions are plot- 
ted in Fig. 3b (green and red line, respectively). At inter- 
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FIG. 3. (a) HAADF-STEM image and ELNES of Ti-L2,3 
and O-K edges across the (6|6) superlattice. (b) Simultane- 
ous dark field profile recorded along the line scan (black line) 
and concentration profiles of Sr (green line) and Pb (red line). 
Horizontal lines mark the interfaces, (c) Energy splitting val- 
ues in PTO and STO layers across the superlattice. Dashed 
lines indicate the values for STO substrate and 12 nm thick 
PTO. (d) High resolution Ti-L2,3 and (e) O-K spectra across 
the PTO-STO interface after background substraction. (f) 
Ti-Ls edges obtained after 10 R-L iterations. 



faces, i.e. units cells labeled 3, 9, 15 and 20, both signals 
clearly cross inbetween neighboring atomic columns, con- 
fining any possible Pb/Sr interdiffusion to ±1 unit cell 
from the interface. 
To relate the observed continuous evolution of the STO 
and PTO spectra across the superlattice to structural 
variations, density functional theory within the local den- 
sity approximation (LDA)^® and charge transfer multli- 
plet calculation^ were performed. Two sets of struc- 
tural models, one for STO and one for PTO, were built 
by artificially varying the c/a ratio from 1 to 1.063, cor- 
responding to the limits of bulk materials. To repro- 
duce the experimental in-plane epitaxial constrain, we 
set a = 6 = 3.905 A. 

Figures 4a and 4b display the relaxed Ti and O atomic 
positions obtained for PTO and STO unit cells with dif- 
ferent c/a ratios, respectively. For PTO, excellent agree- 
ment was found between the experimental and relaxed 
atomic positions of the room temperature PTO unit cell 
(marked with a square]!^. Note that as the c/a ratio is 
reduced, the atomic position of Ti, 0(1) and 0(2) re- 
main strongly displaced from the ideal centrosymmet- 
ric positions, except for the cubic unit cell (c/a = 1). 
Our LDA calculations for STO confirm the paraelectric 



state c/a = 1 as the predicted ground state but, in this 
case, significant Ti off-centering is only observed once 
the tetragonality exceeds ^ 1.02. From then on, the Ti, 
0(1) and 0(2) atoms move further away from the cen- 
trosymmetric positions, but the resulting polar distortion 
is much smaller than that in PTO with the same c/a. 

The electronic density of states was computed for the 
whole relaxed geometrjfSS!. The projected unoccupied 
density of states on Ti atoms was further projected onto 
the Ti d orbitals and decomposed according to the lo- 
cal symmetry €4^, i.e., d^2^ dx2_y2, d^y and d^z + dy^^. 
This non centro-symmetric point group is induced by the 
ferroelectric displacement of the Ti atom along the c axis. 
The resulting values obtained for the LDA CFS (CFS = 
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y2)-Fi{dj.y, d^z+dyz)) for STO and PTO unit 



cells with different c/a ratio are shown in Fig. 4c and 4d, 
respectively. 

These CFS together with the tetragonal distor- 
tion energy differences {Ai—E{dxy)-E{dxz+dyz) and 
A2—E{dz2)-E{dx2_y2)), displayed in figures 4c and 4d, 
and the estimated band widths obtained by LDA for the 
STO and PTO were used as input data for charge trans- 
fer multiplet calculations, yielding a qualitatively good 
match with the experimental measurements (Fig.4e). In 
particular, the difference between the energy splittings of 
PTO and STO (AE(L3)sto-AE(L3)pto) is computed to 
be 0.60 eV, close to the measured value 0.65 eV. More- 
over, the asymmetric broadening in the Ls-Cg band of 
the PTO simulated-spectra with respect to STO spec- 
tra agrees with that experimentally observed in the PTO 
layers (spectra 9 to 12 in Fig.3f). 

We have also performed similar calculations for PTO 
and STO structures with a tetragonal distortion but forc- 
ing a paraelectric state (atomic displacements ^ = 0) , by 
using the Ti centrosymmetric local symmetry D4/1. The 
LDA calculations result in a smaller cubic crystal field 
evolution but much larger tetragonal distortion energy 
that should induce some spectral difference in the EELS 
data. Although the charge transfer multiplet calculations 
of the Ti 1^2,3 edges using ab-initio computed parameters 
are not reliable enough for a quantitative match with 
experiment, it can provide valuable information about 
the spectral differences expected between FE and PE 
state. Figure 4f shows the Ti-L2,3 spectra as obtained by 
the charge transfer multiplet calculations for PE and FE 
tetragonally distorted PTO unit ceU (c/a=1.063). The 
main splitting of the L3 and of the L2 line do not vary 
much, but features appear due to the strong D4?i distor- 
tion parameters. Such features are not observed neither 
in our experimental EELS data nor in previous rep orted 
XAS data of tetragonaly distorted STO or PTd^^lMJ. 
These results thus exclude a possible centrosymmetric 
arrangement of the atoms in the PTO layers, confirm- 
ing their ferroelectric nature, as already revealed by the 
presence of ferroelectric domains. The large cell-to-cell 
variation of the CFS in the PTO layers can thus be at- 
tributed to a variation in tetragonality and polarization, 
showing that the uniform polarization model, frequently 
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FIG. 4. Relaxed atomic displacements as a function of the c/a 
ratio for (a) PTO and (b) STO structural models. CFS values 
for (c) PTO and (d) STO structural models as a function 
of the c/a ratio, (e) Simulated L3 spectra for PTO (red) 
and STO (black), (f) Simulated L3 spectra for a paraelectric 
(blue) and polar (red) tetragonally distorted PTO unit cell 
(c/a=1.063). 



used to (describe monodomain superlattices'^ does not ap- 
ply to our polydomain samples. 

Combining the results of the ab-initio calculations and 
the measurements of the reference PTO and STO spec- 
tra allows us to estimate the local distortions within the 
superlattice layers. The CFS of the central PTO layer in 
the superlattice is 0.3 eV larger than that of the refer- 
ence PTO layer, corresponding to a significantly smaller 
c/a value of about 1.025 which further decreases towards 



the interface. The STO layers show a weaker CFS vari- 
ation, with the four central unit cells all showing split- 
tings around 0.1 cV below that of the cubic STO sub- 
strate. The small tetragonality is consistent with X-ray 
diffraction measurements on pure paraelectric STO films, 
which yield a c-axis lattice parameter of 3.92 A possibly 
signalling a slight off-stoichiometrji^H. The very weak 
distortion in the STO makes it difficult to distinguish 
whether these layers are polar or not but if any polariza- 
tion is present, it is small and any electrostatic coupling 
between the PTO layer^^U jg therefore likely to be weak. 
The structural distortions illustrated in Fig. 3c may arise 
from the inhomogenous strains associated with the do- 
main structure itself. However, we should note that a 
similar reduction of tetragonality at interfaces or surfaces 
of Pb-based ferroelectrics has also been observed in thin 
films without stripe domain^ and hence the precise phys- 
ical origin of these distortions requires further study. 

In this article, we have illustrated that high resolution 
EELS is a powerful alternative technique for studying 
local ferroelectric distortions at the perovskite unit-cell 
scale. The technique is sensitive enough to measure small 
variations in the L2,3-edge down to 50 meV that can be 
semi-quantitatively related to a c/a evolution as small 
as 1%, allowing us to map with ultra high sensitivity 
cell-to-cell variations in local structural distortions and 
revealing an inhomogeneous polarization profile within 
the ferroelectric and paraelectric layers that previously 
has only been predicted theoretically or inferred by indi- 
rect means. Recent XAS measurements reveal that lin- 
ear dichroism is sensitive to the polarization directiorPi', 
suggesting that this information could also be accessible 
through anisotropic EELS. However, achieving high spa- 
tial resolution in anisotropic EELS measurements is still 
a challenging task. 
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